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Abstract 



We have developed a scheme for reducing LIGO suspension thermal noise 
close to violin-mode resonances. The idea is to monitor directly the thermally- 
induced motion of a small portion of (a "point" on) each suspension fiber, 
thereby recording the random forces driving the test-mass motion close to each 
violin-mode frequency. One can then suppress the thermal noise by optimally 
subtracting the recorded fiber motions from the measured motion of the test 
mass, i.e., from the LIGO output. The proposed method is a modification of 
an analogous but more technically difficult scheme by Braginsky, Levin and 
Vyatchanin for reducing broad-band suspension thermal noise. 

The efficiency of our method is limited by the sensitivity of the sensor used 
to monitor the fiber motion. If the sensor has no intrinsic noise (i.e. has un- 
limited sensitivity), then our method allows, in principle, a complete removal 
of violin spikes from the thermal-noise spectrum. We find that in LIGO-II 
interferometers, in order to suppress violin spikes below the shot-noise level, 
the intrinsic noise of the sensor must be less than ~ 2 x 10 _13 cm/v / Hz. This 
sensitivity is two orders of magnitude greater than that of currently available 
sensors. 
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I. INTRODUCTION 



Suspension thermal noise is among the major sources of noise in the Laser Interferometer 
Gravitational Wave Observatory]] (LIGO) |L|. In mature interferometric gravitational wave 
detectors, such as those of "LIGO-II" interferometers, it is predicted to be significant in the 
broad frequency band of 30-110Hz and to dominate other noise sources in the narrow peaks 
around frequencies of the standing wave modes (so-called violin modes) of the fibers on 
which the LIGO test masses are suspended. Fundamentally, suspension thermal noise arises 
from randomly fluctuating stresses of thermal origin in the suspension fibers. Currently, 
there are three general strategies for reducing the suspension thermal noise: 

1. Reduce the mechanical losses in the test-mass suspension, thereby decoupling the 
suspension fiber's motions from the suspension's thermal reservoir and thence reducing the 
fluctuating stresses in the suspension fibers. There is a current vigorous experimental push 
in this direction (see, e.g., ||, |J, JJJ). 

2. Cool the suspension system ||, or some of its components^]. The above two strategies 
will not be discussed in this paper. 

3. Use displacement sensors to monitor the motion of the suspension fibers, thereby 
recording independently the random forces responsible for the Brownian motion of the test 
mass (Langevin forces), and then optimally subtract this recorded force from the LIGO 
readout, thus achieving a partial compensation of the thermal noise. The general idea 
of thermal noise compensation has been widely discussed (See e.g.s Refs. |7[] and ||). A 
concrete version of it for broad-band suspension thermal noise in LIGO was conceived and 
analyzed by Braginsky, Levin, and Vyatchanin (BLV) ||. Strigin and Vyatchanin have 
recently suggested a scheme which in principle allows analogous compensation for internal 
thermal noise f9|. 

A related approach to thermal noise compensation has been developed by Heidmann 
et. al. [0] and Pinard et. al. [11]. They have demonstrated — both experimentally and 



theoretically - - a partial compensation of high-frequency internal thermal noise by im- 
plementing in hardware the subtraction of the Langevin force from the readout through 
a negative feedback loop. The result is an effective "cooling" of the mirror's mechanical 
motions []. 

In BLV, the broad-band 30-110 Hz component was targeted for reduction using a scheme 



1 All results presented in this paper are equally applicable to LIGO's international partners: 
VIRGO, GEO-600, TAMA, etc. 

2 A combination of cooling of the top suspension point, and adjusting location of the laser beam 
on the test mass to neutralize the thermal noise originating at the bottom of the suspension fiber, 
may significantly reduce suspension thermal noise in the broad frequency band Q. 

3 It seems to us that Heidmann et al's [ 10 1 direct "cooling" technique requires a reference mirror 
which is much "quieter" than the mirror used in LIGO, so it is not clear how practical such an ex- 



perimental arrangement is for LIGO. Nonetheless, [ 10 1 is an interesting experimental demonstration 



of the principle of dynamical thermal noise compensation. 
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where the motion of a suspension fiber is monitored along all its length, and an optimized 
average displacement is recorded and subtracted from the LIGO data. This scheme would 
be difficult to implement in practice, due to the fraction-of-a-micron proximity of the optical 
waveguide-based sensor and the fiber. In this paper, we instead propose to track the fiber's 
motion at a single point, which should be significantly easier to implement. This single-point- 
tracking does not allow significant reduction of the broad-band component of the suspension 
thermal noise, but when combined with optimal subtraction it should be quite effective in 
removing violin peaks from the thermal noise spectrum. The possibility of compensation of 
the thermal noise at violin resonances was briefly pointed out by Pinard et al. [ffTJ in the 
context of the "cooling" technique. 

Distinct from reduction of the suspension thermal noise, there are viable strategies to 
filter out narrow thermal violin-mode noise from the LIGO data [|12|]. By contrast with our 
scheme, a procedure of this kind would filter out both the noise and the signal around the 
violin spikes. This works well in most cases, since violin modes dominate other noise in 
very narrow frequency bands, and it is not likely that the signal in these particular bands 
would be significant. However, if the shot noise is reduced by increasing the laser power 
or by narrow-banding the interferometer's response, the broad bases of violin spikes could 
become dominant, and the genuine thermal noise compensation scheme, such as proposed 
here, could well be of use. 

This paper is structured as follows: in Sec. II we give an intuitive introduction to the 
origin of suspension thermal noise and explain qualitatively why our method of thermal 
noise compensation should work. 

In Sec. Ill we formally introduce a new readout variable, x rea dout = ^tcstmass + ct^fiber, 
where x te stmass and ^fiber are displacements of the test mass and of the independently mon- 
itored point of the fiber, and a is a frequency-dependent number which must be chosen in 
an optimal way. We then show how to use the Fluctuation-Dissipation theorem to compute 
the thermal noise in this new readout variable. 

Section IV presents results of our calculations of the optimized parameter a and of the 
reduced thermal noise. We find (see Fig. that on resonance the optimal value of a is 
(l/ir)m/M, where m and M are the masses of the suspension fiber and the test mass, 
respectively. 

In Sec. V we deal with the practical issue of the imperfection of the sensor which is used 
to monitor the suspension fiber. We find that if the sensor is infinitely precise (i.e., if it 
introduces no noise to the measurement), then the violin spike could be removed completely 
from the spectrum; however, in most realistic situations the effectiveness of our method 
will be severely limited by the sensor noise. For example, in order to bring down a violin 
spike to the shot-noise level at the violin frequency, one needs a sensor with the sensitivity 
~ 1.8 x 10~ 13 cm/vTIz [cf. Eq. (|3lD of the text]; this sensitivity is significantly better than 
~ 10 _11 cm, which is the best sensitivity of currently used interferometric sensors of fiber 



motion [16|. Perhaps the situation will change when the next generation of displacement 



sensors comes to life. 
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II. MOTIVATION AND INTUITION 



In this section we discuss intuitively the origin of suspension thermal noise and how to 
compensate it by monitoring devices. 

Suspension thermal noise can be traced to randomly fluctuating stresses which are located 
in each suspension fiber. These stresses make the suspension fiber move in a random way; 
the fiber will make its test mass move randomly as well by pulling it sideways. It is this 
resulting random motion of the test mass that is referred to as the "suspension thermal 
noise" . 

Clearly, the motion of the fiber and the test mass are not independent. The fiber's 
random displacement is ~ M/m greater than that of the test mass, where M and m are the 
masses of the test mass and the fiber respectively. BLV have argued that by monitoring the 
random motion of the fiber one gains information about the random forces acting on the 
test mass; one can then use this information to effectively reduce suspension thermal noise 
by orders of magnitude, at least in principle. Their scheme requires to measure x^ eT , which 
is a horizontal displacement approximately averaged over the fiber, then construct the new 
readout variable 

•^readout ^-testmass Q^fiber; (1) 

where a is an optimization factor of order m/M that must be calculated theoretically. It 
is this new readout variable which has low thermal noise: The gravitational- wave signal is 
virtually unchanged by using this new readout variable, since a is very small (8 x 10~ 6 for 
LIGO). 

The BLV scheme facilitates reduction of suspension thermal noise in the broad-band 
range of frequencies between the pendulum and the first violin peak, and if needed between 
the violin peaks. At the same time, their scheme has the disadvantage of requiring precise 
monitoring of the motion of the suspension fiber along all of its length^; at the moment 
there is, to our knowledge, no active experimental effort in this direction. 

In this paper we propose a modification of the BLV scheme. The new scheme is easier to 
implement experimentally because it requires monitoring the fiber at a single point some- 
where near the middle, rather than over all of its length. The requirement for the sensor 
sensitivity is also significantly less stringent than in the BLV proposal. However, single-point 
monitoring is not effective at removing the thermal noise in the broad-band region. Instead, 
our scheme is designed only to remove the thermal noise around the "violin spikes", which 
pierce the photon shot-noise floor at frequencies above ~ 300Hz for LIGO. 

If we restrict our vision to a narrow band of frequencies close to some particular violin 
resonance, we will find that the motion of the fiber is remarkably simple: it looks like a 
standing wave with an integer number of half-wavelengths fitting the length of the fiber. 
The amplitude of the wave is a thermally fluctuating variable, but the shape does not 
fluctuate much. By sensing the fiber motion at a single point, we can monitor the fluctuating 
amplitude of the violin-resonance standing wave, and thus infer the random force with which 
the fiber pulls sideways on the test mass, at a frequency close to the violin resonance. 



Of course, one must monitor all of the suspension fibers. 
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As an example, consider the dynamics of a suspension fiber and its test mass close to 
the first violin resonance, which corresponds to a violin mode with half-wavelength equal to 
the length of the fiber. The fiber's horizontal displacement is given by 

Xf(z) = Asin(kz) sin(o; v t), (2) 

where A is the (random) amplitude of the standing wave, z is distance from a point on the 
fiber to the top of the fiber, k is the wave vector, such that kl = tt (with I the length of the 



fiber), and oo v = (ir Mgl/m is the angular frequency of the first violin resonance. The 
force with which the fiber pulls sideways on the test mass is given by 

Ft = -Mg = Mgirj sm(u Y t). (3) 

Assume for simplicity that the test mass has no tilt degree of freedom. This is actually 
the case for the currently planned 4-fibers suspension on each LIGO-II test mass. Then the 
laser beam spot measuring the test-mass position is displaced by the same amount as the 
test-mass center of mass: 

Ff TYl 

^testmass = -7^-^ = —Asm{u y t). (4) 

Muj^ ttM 

We neglect the test mass radius compared to the wire length. Equation (d) relates the 
thermal motion of the fiber, which is almost completely determined by fluctuations in A, to 
the thermal motion of the test mass, at a frequency close to oo v . If we choose to monitor the 
fiber at its midpoint, i.e., Xfiber = %f(z = 1/2), then the linear combination 

m 

•^readout ^testmass TT-^fiber (5) 

is expected to have significantly reduced thermal noise at the first violin peak; cf. Eq. ([!]). 

How much can the thermal noise be reduced when this new readout is introduced? What 
is the optimal readout variable for frequencies slightly off the violin resonance? In the next 
subsection we use the Fluctuation-Dissipation theorem to answer these questions. 



III. THE NEW READOUT VARIABLE: ISSUES OF PRINCIPLE 

For simplicity, let us assume that the test mass is hanging on a single fiber, and has 
no tilt degrees of freedom; see Fig. |l[ In LIGO-II, the test mass is suspended on two fiber 
loops, which from the point of view of suspension thermal noise is equivalent to the test 
mass suspended on four independent fibers, and in this configuration the no-tilt assumption 
is realistic. The generalization of our scheme and analysis to a multi-fiber suspension is 
straightforward. 

We assume that the gravitational- wave interferometer monitors the displacement of the 
test mass Q Xtestmass, and that an independent sensor S measures the displacement of the 
suspension fiber x^ er at some point z . 



5 Here it is assumed that the test mass moves as a rigid body, that is, we neglect the internal 
thermal noise due to fluctuations of the test-mass shape. 
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The observer should then record a new readout variable, 



•^readout -^testmass C^fiberj (6) 

where a is a frequency-dependent coefficient which has to be chosen so that the thermal 
noise in x rea dout is minimized. 



The Fluctuation-Dissipation theorem [13 allows one to calculate the spectral density 



S x (f,a) of the thermal noise in x rca dout, where / is the frequency; see Refs. |14| and || 
for the general method of calculation, and see for the first direct application of the 
Fluctuation-Dissipation theorem to computation of suspension thermal noise. First, we 
imagine that a generalized oscillating force F(t) = F cos(W) conjugate to x re adout is applied 
to the fiber + test mass system (here and elsewhere in this paper u> = 2nf). We introduce 
such a force via the interaction Hamiltonian, 

-f^int — — ^(O^readout- (?) 

As discussed in ||, and as is apparent from Eqs. (|B|) and (|7|), applying the generalized 
force F is equivalent to applying two simple "Newtonian" forces simultaneously: one, with a 
magnitude F(t), to the test mass, and the other, with a magnitude aF(t), at a point z = zq 
on the fiber. 

The next step is to calculate the motion of the test mass and the suspension fiber under 
the action of the generalized force, F(t), and to work out the average power, Wdiss, that 
would be dissipated as heat as a result of such motion. The thermal noise in x rea dout is then 
given by 

8k B T W diss 

S x {f,a) = — (8) 



UJ 







where fee is the Boltzmann constant, and T is the temperature of the suspension fiber (cf. 



Eq. (3.10) of [|TJ). 

The last step is to choose a such that S x (f, a) is as small as possible. 

The oscillatory motion of the fiber+test mass under the action of the generalized force 
F is shown in Fig. [I]. The fiber is strongly bent near three points^: at the bottom and 
top attachment points, and at the point on the fiber where we imagine applying aF (the 
location of the independent sensor). It is near these three points the regions in which the 
fiber bends and heat is produced by the fiber+test mass motion. 

In order to compute the dissipated power W^iss, and then to use it in Eq. (|]) to calculate 
the thermal noise, one must specify a model for dissipative losses in the suspension fiber. 
We assume that the time-averaged power dissipated as heat is given by 

Wdiss = / Ctop^top + Cbottom ^bottom + Cmiddle^middle ) (9) 



6 The length over which the fiber is bent is given by A = \J JE/Mg, where J is the geometric 
moment of inertia for the fiber crossection (it equals to 7T(i 4 / 64 for the circular crossection of 
diameter d), and E is the Young modulus of the fiber material. For typical LIGO parameters A is 
a small fraction of a centimeter, much lass than the length of the suspension fiber 
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where top , ^bottom and 6* m iddie are the amplitudes of the oscillating top, bottom and middle 
angles respectively, see Fig. [I], and ( top , Cbottom and Cmiddie are frequency-dependent constants 
which are determined by the dissipation mechanism. If the source of dissipation is distributed 
homogeneously, then 

Ctop = Cbottom = 2Cmiddle = C (10) 

The amplitudes of the oscillating angles 9 top , bottom and ^middle are worked out in Ap- 
pendix A: 







top 



Fa 
M 



k [1 + a cos [k (I — zq)]] — (oc/g) u 2 sin [k (I — z )} 



'bottom 



M 



gk cos (kl) —uj 2 sin (kl) 

k cos (kl) + (a/g) u 2 sin (kz ) 
kg cos (kl) —uj 2 sin (kl) 



and 



e 



middle 



uF 
Mg 



(11) 



(12) 



(13) 



Here M is the mass of the mirror, m is the mass of the suspension fiber, I is the fiber length, 
z is the distance along the fiber from the top attachment point to the sensor, 



u/^Mgl/m, 



(14) 



is the wavenumber of the standing wave induced in the fiber, and z is the distance from 
the suspension top to the sensor. 

The expressions for 9 top and bottom do not take into account damping, and hence they 
diverge when u approaches the violin-resonance angular frequency. We correct for this by 
replacing uj and k by u(l + i/Q„) and k(l + i/Q n ) when u is close to the n-th violin- 
resonance[], and by taking the absolute value of the now complex expressions for 9 top and 
^bottom- Here Q n is the quality factor of the n-th violin mode; it can be shown, using Eqs. @ 
and (0), that 

nMgl 



Qr. 



2C(/n 



(15) 



for the case when the damping is homogeneously distributed in the fiber and is characterized 
by a single parameter C(/)- 

We can then use Eqs. (|TTJ), ( |T2"D and ([13]) to compute the dissipated power, H^diss, from 
Eq. (0). To minimize the readout thermal noise S x (f,a), we must choose the optimal 
ot = «opt which minimizes Wdiss, i-e. dW^ ss /da = 0. 



We can justify this procedure by decomposing the motion of suspension fiber into normal modes, 
and considering how each mode is driven by the generalized force. When the driving frequency 
is close to a proper frequency of some violin mode, we can neglect excitation of other modes. 
Complexifying the frequency in the way described above is the standard procedure for finding the 
response of a damped harmonic oscillator close to its resonance frequency. Here by "close" we 
mean within a few widths 7 n = u) n /Q n of the n-th violin resonance; in our numerical evaluations 
we used the complex replacement within 107 n from n-th resonance. 



7 



IV. THE NEW READOUT VARIABLE: RESULTS 



When we carry through the calculations outlined in the previous section, we find that 
the optimal readout variable which minimizes the thermal noise is 



x 



opt 



"I - "opt^fib 



(16) 



where 



"opt 
A u 

B n 

B\2 
B13 



A n /(B n + B 12 +B 13 ); 

—2gk jgfccos [A; (I — z )\ + uo 2 [cos (kl) sin (kz ) — sin [k (I — Zo)]]\ 
g 2 k 2 cos 2 (kl) + 2g 2 k 2 cos 2 [k (I - z )) 



to 



2a/ 



to 2 sin 2 (kl) - gk sin (2kl) + 2tu 2 sin 2 [k (I - z )] 



lo 2 sm 2 (kz ) — gksm[2k(l — z )] 



(17) 
(18) 



The optimal value a opt is plotted in Fig. |3] for the case when the sensor is monitoring the 
midpoint of the fiber, i.e. z = 1/2. As is clear from the figure, a opt is frequency-dependent, 
and one should keep this in mind while designing the numerical procedure for the data 
analysis. On resonance (where our method is expected to be most efficient), we get 



a 



opt 



m 



ttM 



6.36 x 10" 



m 



0-3 g/ 



15kg N 



(19) 



which is in exact agreement with our more intuitive calculation in Eq. (|5]). 

The thermal noise with and without the optimal monitoring, S x (f, a opt ) and S x (f, a = 0) 
are plotted in Fig. (solid and dashed lines respectively). We see that the optimal monitoring 
removes completely the first and other "odd" violin spikes from the suspension thermal noise 
spectrum. "Even" violin spikes can also be removed if the sensor is positioned off-center of 
the suspension fiber. 

How precisely do we need to choose a, so that significant noise reduction is achieved? 
The spectral density S x (f,a) is a quadratic function of a. We can use the fact that it has 
a minimum at a = a opt to write 



S x (f, a) = S x (f, a opt ) + AS(f) (^^sA 

V "opt / 



(20) 



where AS(f) = S x (f,a = 0) — S x (f,a opt ). When the noise reduction is effective, AS(f) 
S x (f, a = 0). From Eq. (|20|) we see that so long as 



a 



a 



opt 



< 



Sxjjj "opt) 

\S x (f,a = 0Y 



(21) 



the thermal noise reduction will not be seriously compromised. If the above condition is not 



satisfied, then the thermal noise is reduced by a factor 
only on how well the observer is able to tune a. 



(a 



oi, 



opt 



) 2 /a 2 pt which depends 
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V. INFLUENCE OF THE SENSOR NOISE 



So far we have assumed that the sensor, which is independently monitoring the fiber 
displacement, has an infinite precision, i.e. we have assumed that it does not introduce any 
noise of its own into the readout variable. In real life no sensor is perfect. The sensor readout 
is, in general, given by 

■^■sensor ^fiber ^si (^^) 

where n s is the random noise introduced by the sensor. We assume that n s is a Gaussian 
random variable with noise spectral density, N s (f). 
The total readout variable is then 

■^readout -^mirror '^(■^fiber ^s)j (^^) 

and the total readout noise is then 

■W/, a) = S x (f, a) + a 2 N s (f). (24) 

A new optimal value a new of a has to be evaluated, so that dS tota ,\(f, a) /da = 0. Using 
Eqs. (|20| ) and fl24"D , we get after straightforward algebra: 

AS(f) 

aneM) = AS(f) + aUf)N s (ff opt{f) - (25) 
The optimized noise in the total readout variable is then 

StotM, "now) = S x (f, « opt ) + Harm{AS(/), N s (f)a 2 opt (f)}, (26) 

where Harm{A, B} = AB/(A + B) is a harmonic mean of A and B. We use the above 
equation to compute the solid curve in Fig. |[ 

In the neighborhood of the violin resonance, AS(f) ~ S x (f,a = 0). The criterion for 
achieving significant noise reduction is then 

N s (f)^S x (f,a = 0)/a 2 opt (f). (27) 

On resonance fllTl, 



aft n\ 4k B TQm 

S x {f v ,a = = 28 

and a opt = (l/7r)m/M (here n is the order of the violin mode). Then a first requirement for 
the sensor sensitivity^ (that is, the sensitivity for which significant noise reduction can be 
achieved on resonance) in Eq. ( p7|) becomes 



"Essentially, this requirement amounts to what one would expect intuitively: the sensor must be 
able to resolve the thermal motion of the fiber. 
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7 V m / \10 6 J \ fv J V300K7 

(29) 

Most sensors are broad-band, and their sensitivity is essentially /-independent for the 
range of detectable gravitational- wave frequencies. The precision of monitoring specified 
by Eq. (|29|) is achievable by currently available shadow sensors, which can measure with 
precision as high as 10 _9 cm/v / Hz; therefore violin spikes can be reduced somewhat even 
using current technology. 

However, we would ideally like to reduce the thermal violin spike below the shot noise 
level. Shot noise will dominate over all other noise sources except its violin spikes, at the 
violin-spike frequencies. When one represents it by a random displacement of an individual 
test mass, this shot noise for LIGO-II interferometers is [TJ[] 

/ 



Sshotnoisc(Z) ^ 1 x l(T 18 cm/vHz x I ) . (30) 



The sensor sensitivity requirement in Eq. (B7j) then becomes 



y/N s {f v ) < ^/S shotnoise (fv)/a opt ~ 2 x lO^cm/v 7 !^. (31) 



Ageev, Bilenko, and Braginsky |16| have recently built an interferometer which can measure 
a fiber displacement with precision of ~ 10~ n cm/vTIz. They are hoping to improve this 



sensitivity by two orders of magnitude in the near future |L9|. This would, in principle, 
make an interferometric sensor viable for our scheme. 

However, in a narrow-band search the sensor-noise requirement is more severe than 
Eq. (|31~D by a factor ~ 5, because the shot-noise spectral density is lower. 



VI. DISCUSSION AND CONCLUSIONS 

Among possible sources of periodic gravitational waves, the most promising are neutron 
stars with spin periods of 1 — 5 milliseconds |HJ|]. The frequencies of these stars' gravi- 
tational waves lie in the range which is pierced by the thermal violin spikes. Detection of 
these gravitational waves by LIGO will require changing the reflectivity and the location 
of the signal recycling mirror of the interferometer in such a way that the shot noise is 
reduced significantly in a narrow band around the gravitational- wave frequency. If a violin 
spike happens to be nearby, eliminating it might help the gravitational-wave detection and 
measurement. Our method allows one to completely eliminate violin spikes from the ther- 
mal noise spectrum, at least in principle. Figure. § illustrates this by showing a suspension 
thermal noise curve with and without compensation. 



9 This statement is based on current astronomical observations which use electromagnetic waves 
(light, radio waves, x-rays, etc.). Naturally, LIGO might change our notion about what the most 
promising gravitational-wave sources are. 
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In practice, though, we seek a non-trivial experimental design in which each suspension 
fiber is monitored at a single point by an independent sensor. High sensor sensitivity is 
the key for the method to be effective. For LIGO-II we need the noise spectral density 
introduced by each sensor to be less than ~ 2 x 10~ 13 cm/v / Hz; this would allow us to reduce 
the violin spike below the broad-band shot-noise sensitivity curve. In order to "bring" the 
violin spike below the narrow-band shot-noise sensitivity curve, we need even higher sensor 
sensitivity, by a factor of ~ 5. Two types of sensors are currently being used for detection 
of the fiber displacement: the shadow sensors and the interferometric displacement sensors. 
The shadow sensors of ~ 10 _9 cm/v^Hz have been demonstrated, and 10~ 10 cm/v / Hz should 



be possible ||21|| ; but this still is a far call from what we need for effective thermal noise 



compensation. Interferometric sensors, such as those developed at Moscow State University, 
look more promising. Ageev, Bilenko, and Braginsky |16|] have achieved a sensitivity of 
~ 10~ 11 cm/v / Hz for a sensor based on reflecting light from a polished spot on a steel wire. 
The current experiment in Moscow hopes to achieve significantly higher sensitivity by 
attaching a tiny mirror at the center of the fibeiQ. It remains to be seen whether future 
interferometric sensor designs will allow our method to be practical for LIGO. 
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APPENDIX A: MOTION OF THE PERIODICALLY DRIVEN SUSPENSION 

FIBER 

In this appendix we derive the expressions for # top , bottom and ^ddie, Eqs. (|TTD , (|12|). 
and (|I3|) of the main text. These are the amplitudes of the oscillating top, bottom and 
middle bending angles respectively (see Fig. |l|), when a periodic force F = F cos(ut) is 
applied to the mirror, and simultaneously a force aF is applied to the point z = z on the 
fiber at which the displacement sensor makes its measurement. In our derivation we follow 
the spirit of Appendix A in BLV. 

For convenience, we complexify the driving force: 

F = F e iu)t . (Al) 



10 The mirror might introduce extra mass and mechanical friction to the center of the fiber, which 
would increase thermal noise in our readout variable. Techniques developed in this paper can be 
readily used to analyse theoretically thermal noise for such experimental set-up. 
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The equation of motion of the fiber with a force otF^e applied at z = zq is 



d 2 xi 2 d 2 x { aF 



5(z- z )e 



iuit 



(A2) 



dt 2 dz 2 p 

where Xf(z, t) is the horizontal displacement of the fiber at point z and time t, c is the 



velocity of a wave in the fiber, c = J glM/ m, M is the mass of the test mass, m is the mass 
of the fiber, I is the length of the fiber, g is the earth's gravity, and p is the fiber's mass per 
unit length. 

We look for a solution of Eq. (|A2|) in the form Xf(z,t) = u(z)e lu)t . Then for < z < zq 
we have 



and for z$ < z < I we have 



u(z) 



u{z) = Asm (kz) 



B sin (kz) + C cos (kz) 



(A3) 



(A4) 



where A, B, and C are constants to be determined, and k = uj/c is the wavenumber of the 
standing wave induced in the fiber.. The jump conditions at z = z are: 



u(z + e) 
du 
dy 



Applying these jump conditions to Eqs 



= u(z - e), 
du 
dy 

) and 



aFo 
gM' 



(A5) 
(A6) 



we get 



A = B + C cot(kz ), C = j^j- sm(kz ). 



(A7) 



To close this system of equations, we use Newton's second law for the test mass, projected 
onto the horizontal axis [cf. Eq. (35) of BLV]: 



uj 2 u(1) + g 



'du(z) 
dz 



M' 



(A8) 



By solving together Eqs. (|A7|) and (|A^), we get 



A — — 



1 + a cos [k (I — Zq)] — (a/gk) uj 2 sin [k (I — z )] 
gk cos (kl) — uj 2 sin (kl) 



(A9) 



B 



F Q 1 + a[sin(kl) sin(A;,2 ) + (u 2 j ' gk) sin(/cz ) cos(/c/)] 



M 



gk cos(kl) — uj 2 sin(kl) 



(A10) 



C 



F a 
Mgk 



sin(kz ) 



(All) 
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Using Eqs. 
bent; 



dx 

dtop ~ d~z 



and ( |Al3|) - (|ATl ), we can now work out the fiber's three angles of 
k [1 + a cos [k (I — Zq)]) — (a/g) u 2 sin [k (I — zq)] 



Fo 
M 



■ cos (kl) —uo 2 sin (kl) 



(A12) 



e, 



bottom 



dx 
dz 



z=l 



M 



k cos (kl) + (a/g) u 2 sin (kz ) 
kg cos (kl) -co 2 sin (kl) 



(A13) 



and 



- — 

"middle r% 
OZ 



dx 
dz 



z=zo+e 



&F 
gM 



(A14) 



These are the same as Eqs. (^), (0), and (|13|) of the main text. 
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FIGURES 



FIG. 1. Motion of the suspension and the test mass under the action of the generalized force 
conjugate to the readout variable x readout = Xtestmass + ax &hcr . 

FIG. 2. Suspension thermal noise with and without compensation. The dotted line is the usual 
uncompensated suspension thermal noise. The dashed line is for optimal compensation and a 
sensor that is infinitely precise. The solid line is for thermal noise compensation that is partially 
impeded by the sensor noise. For the case in the figure, the sensor is positioned in the middle of 
the suspension fiber; therefore, only violin modes with an even number of nodes get affected by the 
compensation. By shifting the sensor away from the fiber's middle one could compensate thermal 
noise in the "odd" violin peaks as well as the even ones. 

FIG. 3. Optimized value of a, for the case when the sensor is perfect (dashed line) and when 
the sensor is intrinsically noisy (solid line); cf. Eq. (p5|) . 
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